Giant reflection band and anomalous negative transmission in a resonant dielectric 

grating slab: application to a planar cavity 
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The fundamental optical effects that are at basis of giant reflection band and anomalous negative 
transmission in a self-sustained rectangular dielectric grating slab in P polarization and for incidence 
angle not very far from the Brewster's angle of the equivalent slab, are investigated. Notice, that 
the self sustained dielectric grating slab is the simplest system that, due to the Bragg diffraction, 
can show both the former optical effects. A systematic study of its optical response is performed by 
an analytical exact solution of the Maxwell equations for a general incidence geometry. 

At variance of the well known broad reflection bands in high contrast dielectric grating slab in the 
sub- wavelength regime, obtained by the destructive interference between the travelling fundamental 
wave and the first diffracted wave (a generalization of the so called second kind Wood's anomalies), 
the giant reflection band is a subtle effect due to the interplay, as well as among the travelling 
fundamental wave and the first quasi-guided diffracted one, also among the higher in-plane wave- 
vector components of the evanescent /divergent waves. 

To better describe this effect we will compare the optical response of the self-sustained high 
contrast dielectric grating slab with a system composed by an equivalent homogeneous slab with a 
thin rectangular high contrast dielectric grating engraved in one of the two surfaces, usually taken 
as a prototype for the second kind Wood's anomalies generation. 

Finally, the electromagnetic field confinement in a patterned planar cavity, where the mirrors are 
two self-sustained rectangular dielectric grating slabs, is briefly discussed. 

PACS numbers: 41.20. Jb,42.25.Gy,42.25.Fx 



I. INTRODUCTION 

Diffractive phenomena in dielectric materials are at ba- 
sis of many interesting optical effects shown by the light 
propagation in complex systems. The so called Wood's 
anomalies^, the super-radiant effect^, band gaps in 
photonic crystals 5 , anomalous propagation in left handed 
materials 6 and light localization in amorphous photonics 7 
are some of the most interesting properties that take a 
crucial role in the tailoring of the optical devices, allowing 
a full control of the local photon density of states. 

While the resonant anomalies of optical diffraction 
gratings have been observed at the very beginning of 
modern optics, when they were studied by Wood and 
Rayleigth in the Philosophical Magazine^, all the fur- 
ther interesting effects have been obtained recently, due 
to the improvement in the nano-manipulation of the sam- 
ples. 

It is well knowi^ that in the optical response study 
of a self-sustained dielectric grating slab (SSGS) three 
different photon energy zones can be considered: i) the 
so called equivalent slab model (ESM) zone, from zero to 
the threshold energy of the first diffraction wave, where 
only the zero order propagates, ii) the energy zone of 
the second kind Wood's anomalies (SKWA) that ranges 
from the first diffraction threshold till its escape into the 
vacuum (this zone is also called subwavelength energy 
zone), and finally iii) the zone of the first kind Wood's 
anomalies (FKWA) for higher energies. 

Some years ago two of the present authors pointed out 
that in the SKWA energy range, a very broad reflection 



band (giant reflection band: GRB) can be obtained in 
a self-sustained rectangular dielectric grating slab in P 
polarization for an incidence angle close to the Brewster 
angle of the ES M 1Q i n . They showed that this interesting 
property appears in the presence of a complex coupling 
among travelling, guided, evanescent and divergent elec- 
tromagnetic waves. 

Recently, high contrast sub-wavelength self-sustained 
dielectric grating slabs at normal incidence have 
been used in vertical cavity surface emitting lasers 
(VCLSELs), in substitution of the vacuum/cavity Bragg 
reflector, since they show a much broad reflection band 
and a better response as a function of temperature with 
respect to the massive A/4 Bragg reflector a 12 i 13 . 

It is worth remembering that high reflectivity in sub- 
wavelength dielectric grating slabs can be obtained by 
a destructive interference at the outgoing surface (zero 
transmission) in the energy range between the first and 
the second energy threshold of the diffracted waves 
(SKWA)£ i 14 i 15 . However these high reflection bands are 
usually very narrow in energy except when large Fabry- 
Perot oscillations, due to the zero order diffraction wave, 
are present at low photon energies. An example are the 
rather broad and sharp edged reflection bands that can 
be obtained for S polarized ligh t 10 ! 11 . 

On the contrary for P polarization, close to the Brew- 
ster angle of the equivalent slab, the Fabry-Perot os- 
cillation intensities of the zero order wave are strongly 
depressed and giant reflectivity bands can be obtained 
by a synergic effect among the SKWA sharp resonances 
and the electric field components with high in-plane wave 
vectors j 10 ! 11 . 
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In fact, in a dielectric grating slab, at variance of the 
equivalent homogeneous slab, the Fourier components of 
the electric field with high in-plane wave vectors, can also 
be travelling in the direction normal to the surfaces, and 
a complex interplay with surface evanescent /divergent 
waves can be generated. 

In the present work the conditions for which the giant 
reflection bands can be obtained will be discussed for a 
general incidence geometry. Moreover we will show that 
in the energy range where the first diffracted wave propa- 
gates in the vacuum, an anomalous negative transmission 
can be observed for a positive real dielectric tensor. For 
particular values of the incidence angle all the intensity 
can be carried out by the first transmitted order, that 
propagates in the vacuum with a negative wave vector 16 . 

The aim of the present work is twofold. First of all, 
we will discuss the different optical effects that are at 
basis of giant reflection bands formation in rectangular 
dielectric grating slabs for different optical polarizations. 
Second, by increasing the incident photon energy we will 
study the zone where also the first diffracted wave be- 
comes travelling in the vacuum. In this case a negative 
propagation, that leads to the so called super-lenses ef- 
fect, just observed in 2D-photonic crystal slatfi and in 
2D anisotropic waveguides 17 , can be obtained. Finally, 
the property of the electromagnetic confinement between 
two parallel dielectric gratings in optical micro-cavities 
will be also briefly discussed. 

The plan of the paper is the following. 

In section II we describe the theoretical framework 
used to model the optical response of a self-sustained 
dielectric grating slab. 

The different optical contributions due to travelling, 
guided and evanescent /divergent waves, that are at basis 
of GRBs in dielectric grating slabs for P polarization are 
discussed in section III where a systematic study of the 
optical response for a highly symmetric dielectric tensor 
grating will be presented. We will point out that broad 
reflection bands for a general polarization of the incidence 
ray are a robust property of dielectric grating slabs when 
characteristic resonance conditions are fulfilled. More- 
over, we will underline that "giant reflection bands" are 
obtained in P polarization close to the Brewster angle of 
the equivalent slab model. 

In the same section we study the energy zone where the 
first diffracted wave becomes travelling in the vacuum 
and show that, for selected parameter values, negative 
transmission can be achieved in a slab of ID photonic 
crystal^ by Bragg effects. 

In section IV an interesting application of dielectric 
grating slabs in optical micro- cavities will be suggested 
and briefly discussed. 

Conclusions are given in section V. 



II. THEORY 

The Maxwell equations in photonic crystal slabs for S 
and P polarizations are usually solved as an eigenvalue 
problem with respect to the electric and magnetic fields. 
In the present paper we choose to solve the equations as 
a function of the electric field alone, by magnetic field 
elimination, namely: 



V x V x E(r) = -o-D(r) . 



(1) 



We consider a general scattering geometry of an inci- 
dent plane wave on the SSGS, periodic along the x-axis 
and with wires direction along the y-axis, as shown in 
Fig.l. The plane of incidence performs an angle cp Q with 
the Cartesian plane (x,z) and the angle between the wave 
vector, incident from the vacuum, with the z-axis is $ . 

We model the ID rectangular dielectric gratingii, with 
periodicity d and wires width L x along the x axis, of a 
non-magnetic material (fi — 1), with a bulk local dielec- 
tric function 

e(u) = e'(u) +is"(uj) 
and Fourier transformed dielectric tensor 
e G , G ,(oj) = e 5 G , G > + [Ae(w) + ie"{uj)] f x S(G - G') 

(2) 

where f x = L x /d is the filling factor, Ae(uj) = s'(uS) — e Q 
the dielectric contrast, e the vacuum dielectric constant 
and S{G—G') the Fourier transformed geometrical tensor 
(structural factor): 



S(G - G') 



sin [{G - G')L x /2] 
{G-G')L x /2 



The reciprocal lattice vector components are: G = £27r/d 
and G' = £' 2tt /d with I, £' = 0, ±1, ±2, ±N and N 
oo. 

The limit G —> G' defines the so called " equivalent slab 
approximation" of the grating with dielectric function: 



e (d- L x ) +e(uj)L x _ 
d 

= {l-f x )e + f x s{uj) (3) 



In the energy range of low absorbtion (e"(uj) —> 0) and 
s'{(jj) = e' = const, the resonance conditions among dif- 
ferent in plane wave vectors in the grating slab requires: 

i) a large dielectric contrast Ae = e' — e > 10; 

ii) a large filling factor f x = L x /d = 0.5 1.0 that 
determines also the effective dielectric function e and 

iii) comparable slab thickness and wires width 
L z /L x = 1.0 -T- 1.5 that determine the photon Mie's res- 
onances in each rectangular elementary cell. 

Clearly the simultaneous occurrence of the ii) and iii) 
conditions strongly influences the building up of quasi- 
guided waves in the dielectric grating slab, and their in- 
terplay with the travelling fundamental one. 
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Moreover, while the i) and ii) conditions influence the 
strong coupling among the different in-plane wave vec- 
tor electric field Fourier components, generated by the 
grating periodicity (see eq.(2) where the coupling scales 
as the product of the two quantities), the iii) condition 
determines the destructive interference between the two 
propagating waves (travelling and quasi-guided). 

Notice that the dielectric tensor of eq.(2), with row in- 
dex £ and column index £ f , is symmetric and invariant for 
translation along its principal diagonal, being e its value; 
the matrix elements of all the other diagonals {£' = i ± m 
for m = ±1, ±2, ... ± M) are proportional to the dielec- 
tric contrast value Ae divided the distance m from the 
principal diagonal. Moreover, for commensurate values 
between wire width and grating periodicity (L x /d = v/n 
with v and n integer numbers) the matrix elements of the 
diagonals with m = ±n are zero (see eq.2), as well as all 
the other elements in the diagonals with m a multiple of 
n. 

A systematic study is performed by choosing, as duty 
cycle values of the rectangular dielectric grating, the low- 
est Pitagora's fractions p/(p+l) with integer p=l,2,3,4 
values, and optimizing the L z /L x ratio in order to obtain 
the giant reflection band. 

The dielectric tensor simmetry affects the strong cou- 
pling among the incident zero diffraction (G=0) electro- 
magnetic field component and the higher orders (G ^ 0) 
in the grating slab and it can be directly observed in the 
S-polarized optical response. In fact, for a negligible y- 
component of the photon wave vector (q y — )> 0) and in 
the limit L z — >• oo the Maxwell equations reduce to an 
eigenvalues problem: 

Y,M G , G ,(u)v n (G') = k 2 n ^ n (G) (4) 

G' 

where the dynamical matrix has the same symmetry of 
the dielectric tensor of eq.(2) except for the matrix ele- 
ments of the principal diagonal 

uj 2 uj 2 
M GtG ,(u>) = —e GjG , - q 2 x (G)5 GjG , = —e - q 2 x (G) 

C Cr — Cx C 

that are no more translationally invariant. 

In order to deal with finite-dimensional matrices a 
truncation of the diffraction orders is needed. Consid- 
ering, for example, the simple case with cp Q = in S 
polarization, the choice of the N value to ensure the nu- 
merical convergence of the optical response calculation 
can be achieved by two different numerical approaches. 
We can increase the v- value (y —> N) of the square op- 
tical matrix of eq.(4) (of the order (2v + 1) x (2v + 1)). 
In this case, we increase both the eigenstate basis set 
{y>i(G), y>2(G)...(p2i/+i(G)} and the Fourier basis of the 
reciprocal wave vectors {G} of (2v + l)-dimension. We 
define this usual approach as the square matrix method 
(SMM). 

In a second approach, that we define as the n-diagonal 
matrix method (nDMM), we start with a (27V+1) x (27V+ 
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FIG. 1: (color online) Schematic illustration of the system 
under study: a self-sustained dielectric grating slab, periodic 
in the x direction, infinite and homogeneus in the y direction 
and finite in the z one. 



1) dimensional square matrix in eq.(4), with only the 
principal diagonal matrix elements different from zero. 
Then, we increase symmetrically, with respect to the 
principal diagonal, the number of up and down diago- 
nals with non zero matrix elements, till to reach the nu- 
merical convergence. Notice that the 1DMM defines the 
equivalent slab approximation (ESA). 

In the first case the convergence is obtained with a 
Fourier truncated basis set {<pi(G), <P2(G)...(p2v+i(G)}, 
but with the full interaction between the different in 
plane components taken into account. In the second case 
the Fourier basis set dimension is at convergence, but 
with a reduced order of interactions (off-diagonal matrix 
elements). Since these two different criterions give physi- 
cal complementary informations on the behaviour of pho- 
ton propagation in a dielectric grating slab, we will use 
both in the present work. 



III. OPTICAL RESPONSE 
A. Photonic modes dispersion 

Let us consider an incident electric plane wave of en- 
ergy huj and wave vector q on a surface of a rectangular 
SSGS that accomplish the photonic crystals condition 
£fr(cj) « s r b (u}) ^> s'l(uS). Moreover, in the energy zone of 
very low absorption (s b (uj) — 0) we can make the further 
approximation e' b (u) « const. 

The model computation is performed for a P polarized 
incident wave in the resonance conditions regime. All the 
necessary equations for the optical response calculation 
are given in Appendix A. 

We use as bulk refractive index value = = 3.34, 
filling factor f x = L x /d =3/4 and spatial periodicity 
d=300nm, while the plane of scattering is taken nor- 
mal to the wires direction (the angle of scattering plane 
(f = in Fig.l), and the incidence angle is $ = 60° 
rather close to the Brewster angle of the effective di- 
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FIG. 2: (color online) P-polarized dispersion curves k n (uo) 
(solid black line) of a grating with bulk refractive index 
value rib = 3.34, filling factor f x =3/4 and spatial periodicity 
d=300nm. The incidence plane is normal to the wires direc- 
tion ((p = 0), and the incidence angle, r} = 60°. The disper- 
sion curves of the ESA (dot-dashed blu curves), as well as the 
two lowest dispersion curves q z (G) for G=0 and G — —27r/d x 
in vacuum (dashed red curves) are also shown. 



electric function ($ = 72°). The slab thickness L z = 
351.9nm (L z /L x « 3/2) is optimized in order to have a 
strong coupling among in-plane wave vector components. 

In the dispersion curves k n (uo) of the dielectric grat- 
ing, shown in Fig. 2 (solid black line), three different en- 
ergy zones, just presented in the Introduction: i) equiva- 
lent slab model (ESM), ii) second kind Wood's anomalies 
(SKWA), and hi) first kind Wood's anomalies (FKWA), 
can be easily identified in correspondence of different 
treshold energies. 

Notice, that while the zero order curve (n=l) is trav- 
elling for any photon energy, the curves with n > 1 
show different threshold energies (in the limit k z — >• 0). 
For the chosen parameter values and for photon energy 
Tiuo = 1.2eF, a bit greater than the threshold energy 
of the first diffracted mode (huju ~ 1.08eV), the cor- 
responding wavelength, computed for the equivalent slab 
dielectric function e = f x £b + (1 — fx)^o — 8.6167, is 
A = 2tt jk = 351.882nm rather close to the chosen slab 
thickness A « Lj±. 

The dispersion curves in the ESA (dot-dashed curves), 
computed by taking equal zero the off-diagonal elements 
of the block matrices of eq.(A7), are also shown in the 
same picture for different diffraction orders. Notice, 
that the energy thresholds for the first and the second 
diffracted wave G = =pl are very close to that of the 
grating modes (n=2,3); therefore, a rough estimation of 
the thresholds energies for n=2,3 diffracted waves, can 
be obtained by the propagation conditions of the effec- 
tive homogeneous model: 



Moreover, since the dispersion curves k\r>MM have 
been obtained by the diagonal elements of the dynam- 
ical matrix of eq.(ll), the kmMM(G) and k\BMM ( — 2G) 
curves show a crossing that is removed when the off- 
diagonal matrix elements are considered. 

The two lowest dispersion curves q z (G) for G=0 and 
G = — 2i\ j d in the vacuum are also given in the same pic- 
ture (dashed red curves). The G=0 line identifies, among 
the grating modes, the ones that are confined in the pat- 
terned slab while, for each k, the q z (—27r/d x ) curve gives 
the treshold energy above which the -1st order can be 
observed in reflection and/or transmission. In the next 
section we will point out how the grating parameters can 
be choosen to enhance the -1st order transmitted inten- 
sity while depressing the zero order reflection and trans- 
mission. 

Notice that the threshold energy (fioj = 2.21eV) of 
the G = — 2ir/d wave in the vacuum is rather close to 
the threshold (huj n =3 = 1.97 eV) of the second diffrac- 
tion wave in the grating; moreover, while the G=0 wave 
emerges from the sample with the same angle of the in- 
cident one, the G = —2i\jd ray shows a negative prop- 
agation very sensitive to the photon energy, as will be 
further discussed in the next section. 

In Fig. 3a) and 3b) the G-components of the electric 
field eigenvectors {(p n (G)} with n=l,2,3 of eq.(A7) are 
given for an incidence angle $ = 60°, and photon energies 
a bit higher than the first (huj = 1.2 eV) and the second 
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FIG. 3: (color online) Electric field eigenvectors (f n (G) for 
n=l,2,3 computed for incidence angle $ = 60°, and photon 
energy a) Tiu—1.2 eV b) Tluj =2.1 eV 



(Huj = 1.97eV) threshold energies (see Fig. 2). 

It is well known that for normal incidence (# = 0°) the 
amplitude components of the n-th order diffracted wave 
show even (2n+l) or odd (2n) symmetry^ with respect 
to the G — »• —G transformation (obviously this property 
cannot be observed in the present calculation, due to 
the large value of the incidence angle choosen). More- 
over, amplitudes distribution for n=l,2,3 shows maxi- 
mum values in correspondence with t = 0, — re- 
spectively and gives an estimation of the coupling among 
waves with different in-plane wave vectors. The rather 
broad amplitude distribution around the maxima become 
rather sharp when the eigenvector is computed at its own 
threshold energy (see the distribution of n=2 and n=3 in 
Fig. 3a) and 3b) respectively). 

Notice that the SKWA effect in tick dielectric grating 
slabii is a rather more complicated phenomenon with 
respect to that studied in planar waveguides with a thin 
grating engraved on one of its surfaces, a system usually 
considered a prototype for the SKWA explanation 9 . 

This different behaviour is schematically reported in 
Fig. 4a) and 4b), where the propagation of the fun- 
damental (G=0) and the first order diffraction wave 
(G = —2n/d) inside the patterned planar waveguide 
is compared with the propagation of the correspondent 
in-plane components (n=l,2) in a SSGS. Notice, that 
the SKWA of the two systems are qualitatively differ- 
ent; in fact, while on the bottom surface of the grating 
slab (Fig. 4b) the first diffraction wave amplitude (n=2 
and G = — 2ir/d) is partially transmitted by the conver- 
sion on the fundamental wave component (G=0), due to 
the boundary conditions, on the planar waveguide case 
(Fig. 4a) an internal total reflection is observed. More- 
over, while the SKWA in the structure of Fig. 4a) is given 
by the destructive interference between the fundamen- 
tal (G = 0) and the first diffracted wave (G = — 27r/d), 
in the structure of Fig. 4b) many components, with the 
former two in plane wave vectors, interfere; therefore, 
we expect that the two structures show, in the opti- 
cal response spectra, rather the same resonant energies 
(SKWA) but different line-shapes. As an example let us 
consider the case of photon energy Tiuj w 1.242eV and in 
plane wave vector ^(0) = 2.8854 • l(T 4 a.?i. (tf = 60°) 
where two diffracted waves propagate in the grating with 
wave vectors along the z-axis k<i = 3.37838 • 10~ 4 a.u. 
and k\ = 8.64865 • 10~ 4 a.u. (as can be obtained from 
Fig. 2). Since the two diffracted waves have their high- 
est component amplitudes at G-\ = —2i\jd and Go = 
respectively (see Fig. 3a) it is possible to determine the 
propagation angles inside the grating slab that in this 
case are ipi ~ 18°30 / and cp-i « — 67°30 / respectively 
(see Fig. 4b) for lattice wave vector value G\ = 2ir/d = 
11.08354- 10~ 4 a.u. Therefore, the phase shift associated 
with path length difference between incident wave and 
the quasi-guided wave in the grating slab gives the rela- 
tion: tg (<p-i) = ~ m-£- with integer m. A more 
quantitative analysis of the optical response differences 
between the structures of Fig. 4a) and 4b) will be fully 
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FIG. 4: Schematic illustration of a) the geometry for the 
thin grating engraved on an homogeneous slab and b) the 
self-sustained grating slab (SSGS) configuration. 



discussed in the next section. 



B. Reflectivity and transmittivity 

We would like to remind that in a dielectric grating 
supported on an homogeneous slab (see Fig. 4a) the opti- 
cal properties (reflectivity and transmittivity of the fun- 
damental and diffractive waves) are led by two differ- 
ent Wood's anomalies 9 , namely: diffractive anomalies (or 
FKWA) due to the propagation of deflected rays at vac- 
uum sides of the system and wave guiding anomalies (or 
SKWA), due to the interference between incoming light 
and first diffracted wave, guided into the slab (notice, 
that these two effects are given in a decreasing photon 
energy order). 



1. Giant reflection band and negative propagation 

In order to compare the optical properties of the two 
structures of Fig. 4a) and 4b) we choose, for the engraved 
grating of Fig. 4a), the same parameter values of the self- 
sustained one, with a grating thickness (A = L z /10) one 
order of magnitude lower than the SSGS, and the homo- 
geneous slab thickness L = L z — A with a dielectric con- 
stant e s = e. Notice that the former parameter values are 
chosen in order to observe rather the same Fabry-Perot 
and SKWA effects in the two structures. 

Moreover, while in the former structure (Fig. 4a) the 
source of the diffracted waves and of the quasi-guided 
waves is given essentially to the periodic grating and 
the homogeneous slab respectively, in the latter struc- 
ture (Fig. 4b) the two effects are completely mixed in all 
the slab volume. 

In Fig. 5a) and 5b) the optical reflection spectra, for 
an incidence angle $ = 60° and in the SKWA energy 
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FIG. 5: Reflected intensity, as a function of energy, for a 
self-sustained grating of thickness L z — 351.9nra (solid red 
curve) and for a thin grating (L z /10) engraved on a thick slab 
(9L Z /10) of dielectric function e (dot-dashed green curve): 
a) S polarization and b) P polarization for incidence angles 
(p = and $ = 60°. c) reflection phase for the P polarization 
case, for the engraved (solid red curve) and the self-sustained 
grating (dot-dashed blue curve) structure. 



range (see Fig. 2), for both the structures of Fig. 4a) and 
4b) are shown for S and P polarization respectively. In 
Fig. 5a) we observe that the optical spectra for S polar- 
ization show Fabry-Perot oscillations, for photon energies 
lower than the first treshold, very similar for the two dif- 
ferent structures, and the SKWAs at rather the same 
energies, has hypothesized beforeii. At variance, for P 
polarization (Fig. 5b) narrow and broad reflection bands 
are observed for the engraved and the self-sustained grat- 
ing structure respectively. In this case, since the angle 
of incidence is very close to the Brewster angle of the 



equivalent slab model, the Fabry-Perot oscillation ampli- 
tudes are rather small for the engraved system while are 
completely suppressed for the SSGS and they show very 
different reflection spectra. 

For the SSGS any possible anomalous behaviour, in 
the band where the reflectivity intensity remains equal to 
unity, can be identified only in the phase of the reflected 
wave. To this end the reflection phase versus energy, for 
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FIG. 6: (color online) P-polarized reflected (solid red line) and 
transmitted (dashed blu line) intensity for ip Q = 0°,$ = 60°, 
for a) the G=0 wave in the energy range from zero to the 
second energy threshold and, in the FKWA range, b) the G=0 
wave, c) the first diffracted one. In a) the reflection curve, 
Resm, of the ESM (dotted orange curve) and the reflection 
peaks Rn=i of the SKWAs (dot-dashed green line) are also 
shown. In c) the transmittivity computed considering only 
the G = 0, dz27r/d components in the dielectric tensor is shown 
(dot-dashed green line). 
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the P polarization case, is shown in Fig. 5c) for both 
the engraved and the self-sustained grating structure. It 
can be seen that, as expected, the reflection phase for the 
engraved structure has a change equal to 2tt at the reflec- 
tion peak position. The changes in the reflection phase 
of the self-sustained grating structure are a signature of 
similar resonances underlying the giant reflection band. 

Now we consider the optical response for P polariza- 
tion of the SSGS in the same energy range of Fig. 2. The 
optical spectra for reflected, transmitted and deflected 
rays are shown in Fig. 6. In Fig. 6a) is shown the optical 
response in the energy range from 0.5 eV to the energy 
where the first diffraction wave (G=— 2it/d) escapes in to 
the vacuum (see dispersion curves of Fig. 2). The three re- 
flection curves correspond to three different basis trunca- 
tions in the reciprocal lattice wave vector space, namely: 
i) for G=0, the reflection curve, Resm-, of the ESM is 
obtained, ii) for G = 0, ±2ir/d the reflection peaks Rn=i 
of the SKWAs are shown and finally iii) at convergence 
(G = i 2ir/d for £ = 0, ±1, ±2, ±N for TV « 10), where 
also surface waves are present and the giant reflection 
band appears. The transmittivity is also shown in the 
same picture in order to check the computation accuracy 
(R+T=l). 

In the case i) we observe Fabry-Perot oscillations of the 
equivalent slab model very similar to that of the engraved 
system shown in Fig. 5b), for photon energies lower than 
the first treshold. In the ii) case four well shaped SKWAs 
appear, distributed in all the energy range between first 
and second threshold energies, rather similar to those ob- 
served in the engraved dielectric slab (see Fig. 5b). There- 
fore, at this degree of approximation, the similarity be- 
tween the two different structures of Fig. 4a) and 4b) 
is recovered, with the difference that, at variance of the 
engraved grating structure, in the self-sustained one the 
SKWA peaks are very sharp in correspondence of the en- 
ergy thresholds (n=2 Tiu ~ 1.25eV, n=3 Tiuj ~ 1.97 eV), 
while rapidly broadens for energy far from the thresholds, 
due to the interference of many in-plane components (see 
Fig. 3a and 3b). 

Finally, when the numerical convergence is reached iii) , 
a giant reflection band is obtained close to the Brewster 
angle value with the crucial contribution of higher or- 
der diffracted wavesii. For photon energy Tiuj « 2.2eV^ 
rather close to the treshold of the second diffracted wave 
and to the escape into the vacuum of the first diffracted 
wave, the synergic effect between SKWAs and higher or- 
der diffracted waves, that is at basis of the giant reflection 
band, deteriorates and the reflection band disappears. 

The mechanism that undergo giant reflection bands in 
P polarization, near the Brewster angle, is then rather 
different from that observed for broad reflection bands in 
S polarization (or in P polarization for an incidence angle 
rather far from the Brewster angle). In the S polarization 
case, from the first to the second energy threshold, well 
shaped reflection bands originate due to the interference 
between the large Fabry-Perot oscillations, in the slab 
resonance condition, and the quasi-guided first diffracted 
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FIG. 7: (color online) Output angle i9 ou t of the -1st diffracted 
wave as a function of energy. 



wave (SKWA); higher order waves seem to give a minor 
contribution to the reflection effect (see ref. [10] and 
Fig. 5a). 

Now we move to the energy zone of the first kind Wood 
anomalies and show, in Fig. 6b) and 6c), the reflectivity 
and transmittivity intensities of the fundamental mode 
G=0 and the -1st diffracted one that propagates in the 
vacuum. 

Also in this range of photon energies is easy to verify 
that the sum of all the propagating reflection and trans- 
mission waves equals 1. 

It is worth noting that, in the energy range close to 
the first diffracted wave escape into the vacuum, the de- 
flected ray shows the highest intensity with transmission 
close to the unity (~ 94% for Tiuj = 2.3eV), while the 
zero order reflected and transmitted waves are drasti- 
cally suppressed. Moreover the -1st diffracted wave, with 
Q// = (Qx — G)i propagates in a left direction with an 
energy dependent output angle $ out as shown in Fig. 7 
for all the energy range of the optical spectrum. For 
example, at photon energy Tiuj = 2.3eV the deflection 
angle is ~ 20°, therefore an incident wave will be 
deflected with an angle as large as 50°, as shown in the 
inset of Fig. 7, and intensity greater than 90% as shown in 
Fig. 6c). The transmittivity computed considering only 
the G = 0,±27r/d components in the dielectric tensor, 
shown with a dot-dashed green line in Fig. 6c) ,has a 
rather different behaviour respect to the full calculation. 

In conclusion, we would like to underline that the nega- 
tive propagation is produced by the Bragg effects in a non 
magnetic material with positive bulk dielectric function 6 , 
a completely different effect with respect to that hypoth- 
esized by Veselago many years agoi£, due to the nega- 
tive dielectric and magnetic susceptibilities of the sam- 
ple. Moreover, because this effect is not present in the 
engraved structure under the same scattering conditions, 
we can conclude that also for this property the higher 
order diffracted waves are crucially important. 
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FIG. 8: Reflection intensity as a function of photon energy and in plane wave vector for three different grating with filling 
factors values a) f x = 2/3, L z = 320nra, b)f x - 3/4, L z = 350nra, c)f x - 4/5, L z = 380nra 



2. General scattering geometry 

In order to go a bit deeper in studying the giant reflec- 
tion band and the negative propagation present in two 
different and characteristic energy zones close to the first 
and the second Wood's anomalies, let us consider the 
scalability properties of the SSGS. 

As discussed in ref.[ll] the giant reflection band is scal- 
able in energy if the filling factor f x = L x /d and the ratio 
L z /L x are taken constant. For the present parameter val- 
ues and f x = 3/4, L z scales according to the wavelength 
of the equivalent slab model, computed at the first pho- 
ton threshold energy (L z ~ A) of the grating as shown in 
Tab.l. 



TABLE I: Grating parameters 



uJo(eV) 


Auj(eV) 


L x (nm) 


d(nm) 


L z (nm) 


A(nra) 


1.818 


1.273 


150 


200 


234.6 


234.588 


1.212 


0.848 


225 


300 


351.9 


351.882 


0.909 


0.636 


300 


400 


469.2 


469.176 


0.454 


0.318 


600 


800 


938.4 


938.352 



While the filling factor determines the effective dielec- 
tric function value of the grating, the ratio L z /L x is 
connected with the Mie energy states of the elementary 
nano-particle of the system; therefore, the two basic ra- 
tios L x /d and L z /L x completely determine the building 
up of the quasi-guided waves in the former dielectric grat- 
ing slab. Notice that the giant reflection band can be as 
large as the whole visible energy range, and moreover, 



it scales from infrared to ultraviolet energies according 
to the simple relationship: Auo/uo « L x /d ~ 3/4 and 
L z /L x w 3/2, where the band width Au and the lower 
energy band edge uo are calculated for reflectivity values 
1/2. Notice that for R> 98% the band ratio is AA/A > 
40% where A is calculated at the band center. 

The former properties are not restricted to the parame- 
ter values considered in Tab.l. In fact, in Fig. 8a)-8c) we 
plot the reflection intensity as a function of photon energy 
and in plane wave vector for three different grating with 
filling factors values L x /d = p/{p + 1) where p=2,3,4; 
the ratio 1/2 is not taken into account since it does not 
give a real giant reflection band, due to the small value 
of the effective dielectric constant that cancels the strong 
coupling condition among different in plane wave vectors 
components. The grating periodicity is taken constant 
(d=300nm), the wire width scale as: L x = pd/(p + 1) 
and the slab thickness (L z /L x « 3/2) is optimized in or- 
der to obtain a reflection zone as large as that observed 
for the sample with L x /d = 3/4. In fact rather similar 
zones of high reflectivity (in red in the plots), for the three 
different filling factor values, are present and in energy 
cover the range from the first to the second diffraction 
thresholds. 

We remind that in all the three different cases the di- 
electric matrix of the grating is symmetrical and trans- 
lational invariant with respect to its principal diagonal, 
and indeed has zero value for the up and down diago- 
nal matrix elements far from the principal one with pe- 
riodicities m, where m is the denominator of the filling 
factor: L x /d = n/m (see eq.2). Notice that the simme- 
try strongly reduces the number of non zero matrix ele- 
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FIG. 9: (color online) Reflection (R) and transmission (T) 
intensities of the two electromagnetic waves travelling in the 
vacuum: G=0 and G — —27r/d in the energy range 2.1-1-3.2 eV 
in the second Brillouin zone. The same color scale of Fig. 8 is 

used. 



merits of the block dielectric tensor (27V + 1) x (27V + 1) 
to 27V(m — l)/m reducing the interaction among differ- 
ent in-plane wave components. Since the giant reflection 
band in P polarization is a robust property of the dielec- 
tric grating slab under resonance conditions, it can be 
observed for filling factor values in all the range 0.5 
1.0. 

The folding of dispersion curves at Tiuu = 2.06 eV pho- 
ton energy, at the boundary of the first Brillouin zone, is 
clearly shown in Fig. 8; this energy is close to the escape 
energy into the vacuum of the first diffracted wave (for 
Q// = (Qx — (see Fig. 2). Moreover, for photon en- 
ergy of the fundamental wave (G=0) greater than 2.1eV 
a very small reflection intensity is observed for all the 
three dielectric grating slabs, as observed before in the 
sample p/(p + 1) = 3/4 at $ = 60° incidence angle (see 
Fig.6b-6c). 

In particular in Fig. 9 reflection and transmission inten- 
sities are shown for the two electromagnetic waves travel- 
ling in the vacuum (namely: G=0 and G = — 2ir/d) in the 
energy range 2.1 -=- 3.1 eV in the folded second Brillouin 
zone. 

We observe high intensity value only for the -1st trans- 
mission wave, therefore the optical behaviour of the di- 
electric grating slab is similar to a negative propaga- 
tion for a rather broad energy range and incidence an- 
gles. In fact, in the second Brillouin zone, the intensity 
of transmitted and reflected optical waves are negligi- 
ble small except for the transmission of deflected ray at 



Q// = (Qx — G)i, and this behaviour is in agreement with 
the optical spectra of Fig. 6c). 

In Fig. 10a) and 10b) the reflectivity maps for TE and 
TM polarization respectively are computed for a p/(p + 
1) =3/4 grating slab and angle of incidence $ = 60°, as a 
function of < ip < tt/2 angle (see Fig.l); the reflection 
spectra for ip Q = are also shown for sake of comparison. 
In the central maps are given the reflection intensities 
computed with only the first three waves (N=l) in the 
Rayliegh expansion of the electric field, while the right 
maps are at convergence. 

The TE polarization reflectivity shows a broad zone of 
high reflectivity values and the optical response is rather 
insensitive to the contribution of higher order evanescent 
waves, while for TM polarization where the role of the 
evanescent waves is crucially important the broad reflec- 
tion band still remains very sensitive to the incidence 
angle values. 




90 cp 90 



r 22 

- 2.0 

- 1.8 






' 9 


if 




" 1 ' 4 | 
-1,1 

- 1.0 






V" 

j- 0.6 







Reflectivity o q, 90 o q, 90 

FIG. 10: (color online) Reflectivity maps for a) TE and b) TM 
polarization for a f x = 3/4 grating slab, angle of incidence $ = 
60°, as a function of the incidence plane angle ip (see Fig.l). 
The blu curves on the left, corresponding to the reflection 
spectra for ip Q = 0, as well as the central maps are obtained 
for N=l; the right maps are at convergence. The same color 
scale of Fig. 8 is used. 
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IV. PATTERNED MIRROR CAVITIES 



be written as: 



The giant reflection band effect discussed in the previ- 
ous sections suggest the use of sub-wavelenght gratings 
as very efficient reflectors in vertical microcavities since 
they allow the use of a single patterned layer, contrary 
to multilayer Bragg mirrors, and give wider wavelength 
ranges of high reflection. 



Recently, a self-sustained rectangular dielectric grat- 
ing was used at normal incidence angle in vertical cavity 
surface emitting laser o 12 i 13 . 



The interest for the former system with respect to the 
massive Bragg reflectors is two folds: i) in energy a broad 
stop band is obtained ii) it should show a better response 
as a function of the temperature variation. 



In Fig. 11 the design of a vertical cavity is presented 
where the electromagnetic field confinement is achieved 
by the internal reflection of two patterned mirrors at dis- 
tance L. The case of mirrors with the same periodicity of 
the lateral pattern and the same slab thickness d\ = d2 is 
considered. If we take the in-plane wave vector along the 
x-axis (cp = in Fig.l) the Maxwell equations (eq.(A7) 
for cp = 0) separate in S (eq.(4) and P polarization. 



The in-plane electric field components in the cavity, in 
the grating slab and in the region to its right side, can 



El 


B L 




< 



B, 



FIG. 11: Schematic design of a patterned cavity 



< z < L/2 



iqx (G)x 



A R {G)e 



(G)( z -f)_ 



+B R (G)e 



(5a) 



L/2<z< d 2 +L/2 
E // (r;co)=J2 



A 2 (n)e ik " {z -^- 



+B 2 (n)e~ 

d 2 + L/2 < z < 00 
E,,{r-uj) = Y,e tqAG)x 



ik n (,2 cy ) 



iq x (G)x 



(5b) 



A 3 (G)e zqAG)(z - d2 -2 ) + 
+ B z {G)e iqAG){z - d2 -^ ) 



(5c) 



where the k n 2N+1 square roots of the eigenvalue prob- 
lem, and the {cp n (G)} correspondent 2N+1 eigenfunc- 
tions with 2N+1 G-components of the Maxwell equation 
(eq.A7), are computed for a chosen q x (0) value in the 
first Brillouin zone (-7r/d, n/d) . 

A phase matrix: 

/ e iq z (-N2n/d)L Q \ 



X > (L) 







e «z*(o)L o 
\ e iq ^ N27r ^ d ^ L J 

match the fields in the left and right side of the cavity: 



Ar= X > (L)A l B r = X < (L)B l 



(6) 



where the X < (^) matrix elements are obtained by the 
substitution: e iqz ^ L e~ iqz ^ L . 

By considering the S polarization, and imposing the 
Maxwell boundary conditions at z = L/2 and z = ^2 + 
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L/2 we obtain the amplitudes in the transfer matrix form: 



{ A R (G) = Z[t>(G,n)]- 1 cp n (G) 

n 

{A 2 (n) + r>(G,n)B 2 (n)} 
B R (G) = E[t>(G,n)}- 1 lfin (G) 

{r> n (G,n)A 2 (n) + B 2 (n)} 
{ A 3 (G) = Y,lt 2 > (G,n)}- 1 MG) 

n 

{A 2 (n)e ik " d2 + r^(G,n)B 2 (n)e 



(7a) 



B 3 (G) = J2[t>(G,', 



Vn(G) 



{r$(G,n) e ife " d M 2 (n) 



where: 



r>(G,n) 



q z (G) - kn 
q z (G) + k n 



t>(G,n) 



ik ^B 2 (n)} 



2g z (G) 

Qz(G) + k n 



(7b) 



(8) 



are the interface reflection and transmission amplitudes 
of the right grating slab. 

It is however well known that, as discussed in ref. 
[14,19,20], the inevitable existence of evanescent solutions 
in dielectric gratings makes the transfer-matrix calcula- 
tion numerically unstable very quickly, so that the scat- 
tering matrix formalism must be used. 

In the scattering matrix form the equations 7a) and 
7b) that now read: 



A R (G) = j:[t>(G,n)}- 1 l p n (G) 

n 

{A 2 (n)+r>(G,n)B 2 (n)} 
S 3 (G) = E[* 2 > (G,n)]-Vn(G) 

n 

{r>(G, n) e ik ™ d2 A 2 (n) + e~ ik " d2 B 2 (n)} 



(9a) 



A 3 (G)=J2[t>(G, 1 



fn(G) 



{A 2 (n)e ik ~ d2 + r% (G, n)B 2 (n)e- ik ~ d2 } 
B R (G) = Z[t>(G,n)}- 1 i pn(G) 

n 

{r>(G,n)A 2 (n) + B 2 (n)} 



(9b) 



give the input and output field amplitudes of the right 
grating slab. 

The forward optical response of the dielectric grating 
slab, for a given incident in-plane wave vector in the ex- 
tended Brillouin zone Q// = [q x (0) + G']i = q x (G')i, 
defined by the assumption: 

F) Ar(G) -+ 5 GtG ,, B R (G) -+ r 2 (G,G'), A 3 (G) -+ 
t 2 (G,G f ) and B 3 (G) 0.0 

makes the system of eq. 9a) an heterogeneous algebraic 
system that can be solved with respect to the values of 
the internal electric field amplitudes {A 2 (n), B 2 (n)}. 



Then the system of eq.9b) gives the matrices of forward 
reflection r 2 (G, G') and transmission t 2 (G, G') whose di- 
mension is (27V + 1) x (27V + 1). 

Due to the optical simmetry of the grating, the back- 
ward optical response, defined as: 

B) A R {G) -+ 0.0, Br(G) -+ t 2 (G), A 3 (G) -+ r 2 (G) 
and B 3 (G) -> 5 G ,G' 

give the same values of reflection r 2 (G, G') and trans- 
mission t 2 (G, G') amplitudes of the forward one. 

We can then write the scattering matrix of the system 
as a 2 x 2 block matrix in {G, G'} in the form: 




S 2 (u) 



Afl \ f t 2 (cj) r 2 (u) \ I A R \ 
B 3 J \t 2 (uj) t 2 (uj) J \ B 3 J 

(10) 

An analogous relation holds for the left grating slab: 



Al I = S (lu)( A °\ = ( ri ^ 
Bo i \Bl i I rxH ti(w) 



Moreover, eq. (6) now is: 





X > (L) A, 
X <(L) J 1 B Z 



(11) 



(12) 



where the field amplitudes are given in the {G, G'} re- 
ciprocal space, and X^{L) (f° r P = < •> >) are diagonal 
matrices. 

Moreover, from eqs. (11), (12) and (13) we obtain: 



A L = G>(L) [ti(w)A + r 1 ( w ) X > (f)t 2 ( W )B 3 ] 
B L = G<(L) [r 12 (w)ti(w)A + x>(L)t 2 (w)B 3 ] 



where the tensors G(L) are: 



(13) 



G>(L) = [i-nHnsH]" 1 

G < (i) = [i-r 12 HnH]- 1 (14) 

with r 12 (uj) = x > (L)r 2 (to)x > (L). 

The total scattering matrix of the cavity is then: 




t 2 (w) X > (L) 
ti(w) 

( r 2 (w) N 
I nH , 




(15) 



Notice, that the tensor G(L) embodies the poles of the 
resonant matrix, and, since the poles of the matrix are 
the zeros of its inverse, these are the eigen-energies of 
the electromagnetic field confined between the two grat- 
ing slabs in the cavity. Moreover, the reduced scattering 
matrix embodies also the optical properties of the iso- 
lated grating slabs, as shown in the second term of the 
right side of eq.15). 
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V. CONCLUSION 

The physical effects that are at basis of giant reflection 
band and negative light propagation in a self-sustained 
rectangular dielectric grating slab are clarified by com- 
parison with an analogous optical system obtained by 
an homogeneous planar dielectric waveguide with a thin 
grating engraved on one of the surfaces^. Notice, that the 
latter system is usually suggested in the literature^ as a 
prototype for explaining the SKWA in a dielectric grat- 
ing. The role of the strong coupling among electromag- 
netic waves with different in plane wave vectors on the 
SKWAs and of the symmetry of the dielectric grating ten- 
sor are also discussed. The different physical effects that 
are at basis of broad reflection bands in S and P polariza- 
tion are clarified. We pointed out that in a self-sustained 
grating slab, under resonant conditions, giant reflection 
bands and negative transmission can be obtained in P po- 
larization, due to the interplay among travelling, guided 
and evanescent /divergent waves not very far from the 
Brewster angle. The optical behaviour for a general scat- 
tering geometry and mixed polarization is also briefly dis- 
cussed. Finally, the attitude on the confinement of the 
electromagnetic field on laterally patterned planar cavi- 
ties is also discussed for S polarization. 



FIG. 12: S-polarized reflection spectra for a laterally pat- 
terned planar cavity obtained with two parallel dielectric grat- 
ing slabs of filling factor f x = 3/4, elementary cell ratio 
L z /L x — 3/2 and periodicity d=300nm for incidence angle 
a) <& = 55°, b) $ = 60°, c) = 65°. 



The S-polarization reflection spectra of a laterally pat- 
terned planar cavity obtained with two parallel dielectric 
grating slabs at distance L, filling factor f x = L x /d = 
3/4, elementary cell ratio L z /L x =3/2, and periodicity 
d=300nm are shown in Figs. 12 a)-c) for three different 
incidence angles $ = 55°; 60°; 65°. The cavity thickness 
(L = 826.34 nm ) is chosen for obtaining a A/2-cavity in 
a guided geometry. The reflection spectra of a single pat- 
terned dielectric grating slab are also shown in the same 
pictures; it is interesting to note that the reflection bands 
of this kind of cavity is improved with respect to the re- 
flection band of a single dielectric grating slab. Moreover, 
the cavity peak shifts in energy linearly as a function of 
the incidence angle, due to the small angle range chosen, 
and the half-width at half height is V « 5meV. 

In conclusion, these results underline that self- 
sustained dielectric gratings are well suited for obtain- 
ing patterned planar cavities in guided configuration; a 
more complete analysis will be performed in a subsequent 
paper. 



Appendix A: Optical response of a self-sustained 
dielectric grating 



A schematic illustration of a self sustained dielectric 
grating slab with general incidence plane wave of Tiuj en- 
ergy is given in Fig.l. 

The electric field components (a = x, y, z) in the grat- 
ing region, < z < L z , expanded in plane waves, are: 



E a (p, z\ u>) = E <* (q// + G l k; uj) e i(q // 



(q//+<3i)-p Jkz 



G 



(Al) 



where p is the in-plane vector (p = (x, y,0)). 

The wave vector of the incident electric field 
q = (q//(0),(fe(0)) has q /f = q x i + q y j = 

cos (p Q i + sin cp Q j and z-component q z (0) = 

~ Qo w ^h q Q = ^y^sin^o, while the in 
plane wave vectors of reflected, transmitted and de- 
flected fields are (G) = (q x + G)i + q y j = q x (G)i + q y j 

and q z (G) = \J ^r£ — Qx (G) ~ Qy- Therefore, the polar- 
ization plane of the deflected i — th wave component 
performs a solid angle (pi = arctg ^ — cos^^^/d ) w ^ n 
the (x,z)-plane. 

Moreover, the unit vectors of the incident electric field 
for TE and TM polarizations, defined by the transversal- 
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ity conditions, are: 

Ste = (- sin cp , cos <p G , 0) 

£tm = (cos cp cos $ , sin cos # , sin $ ) . 



(A2a) 



(A2b) 



The Fourier transformed Maxwell equations in mixed 
coordinates, are: 



92 2 «y 



Bx(G,z)+fe(G)g w f^(G,z) + 

. , n 9^M u, 2 

( G ) = —# D * ( G ' z ) 



03 



2 Qx 



dz 

E y (G,z)+q x (G) q y E x (G,z) 
dE z (G, z) _ cg 2 



(A3a) 



-iq y - 



dz 



D y (G,z) 



,2 



(A3b) 



+ [^( G ) + ^]£ z (G,z) = ^ (G,z) 

(A3c) 

and, given the dielectric displacement field in the form: 
D a (G,z) =Y,£G,G'E a (G',z) = e ikz J2zG,G>E a (G',k) 

G' G> 

(A4) 



they reduce to: 



E 



U 2 X 



E x (G',k)- 



+q x (G) q y E y (G, fc) = fc 2 ^ (G, - fc^ (G) E z (G, fc) 

(A5a) 



E 



G' L 



-T £G,G' - <?x (G) <5g,c 



E v (G',k)- 



+q x (G) ^ (G, fc) = fc 2 ^ (G, fc) - fcg^ (G, k) 

(A5b) 



E 



r w 2 



£g,G' - {ql (G) + g 2 ) 6 GtG < 



E z (G',k) 



= -k [q x (G) E x (G, fc) + q y E v (G, fc)] 

(A5c) 

The former system can be solved with respect to the 
z-component of the electromagnetic field: 

E z (G,fc) = 

= -fc 2 M^ G , {q x (GO ^ (G', fc) + q y E y (G', k)\A6) 

G' 

2 

where the matrix of elements Mq,G' — £g,G' — 
(g 2 (G') + g 2 ) Sg,G' is rea l an d symmetric in the pho- 
tonic crystal limit [s"{oS) — » 0). The substitution of 



the z-component of the electromagnetic field in the first 
two equations of the system of eq.(A5) gives the Maxwell 
equations as a generalized eigenvalue problem: 



A A 

A T A 

^xy ^yy 



k 



2 ( ^xy 



fx 
Vy 



(A7) 

where each block of the matrix, of dimension 
(2N+l)x(2N+l), real and symmetric in the photonic 
crystal approximation, is given by: 

A XX (G : G') = ^eg,g' — Qy$G,G' 

A xy (G,G f ) = q x (G f )q y 5 GjG , 

A yy (G, G') = £e G , G , - <£(G')5 G , G , 

B XX (G, G') = q x (G)Mc] G ,q x (G') + 5 G , G / 1 } 

B xy (G,G') = q x (G)M^ G ,q y 

B yy (G, G') = q y M^ G ,q y + 5 GiG , . 

We solve the generalized eigenvalue problem by the 
method of refs.[9,10] to which it reduces when q y = 0. 
We first have to diagonalize the matrix 



A A 

A T A 

■^xy ^yy 



by solving: 



where 



(A9) 



(AfO) 



Then, trough the unitary transformation U that diago- 
nalizes A we define the new quantities ^ n = U T £ n and 
B = U T BU and transform the generalized eigenvalue 
problem into a canonical one. 

In fact with the unitary transformation we obtain: 

A£n = & 2 B £ n 

where, A is the diagonal matrix of the complex eigenval- 
ues. 

Finally, defining 



and 



L = a i/2 c; 



B = A 1 /2ba 1 /2 



we obtain the canonical eigenvalue equation: 

TJ in = B i n . 



(All) 



Notice that matrix A and the renormalized matrix B 
are symmetrical and real in the photonic crystal limit, 
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while they are not Hermitian for a complex dielectric 
bulk value. For e" (uo) ^ the eigenvalues become com- 
plex and the eigenvectors are no more orthogonal. In 
this case, the eigenvalue problem, eq. (A7), must be 
solved togheter with its Hermitian conjugate and the bi- 
orthogonality relations between eigenvectors correspond- 
ing to complex conjugate eigenvalues must be used. 
Finally, imposing the continuity of the in-plane elec- 



tric and magnetic fields at z=0 and z = L z interfaces 
we compute the optical response of the dielectric grating 
slabii. 

In the particular case of cp Q = the off-diagonal blocks 
of both A and B are zero and ~B yy reduces to the identity 
matrix. Consequently the system in eq.(A7) separates in 
two eigenvalue problems for S and P polarization. 
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